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We obtain, using semi-analytical transfer operator techniques, the Edwards thermodynamics of a 
one-dimensional model of blocks connected by harmonic springs and subjected to dry friction. The 
theory is able to reproduce the linear divergence of the correlation length as a function of energy 
density observed in direct numerical simulations of the model under tapping dynamics. We further 
characterize analytically this divergence using a Gaussian approximation for the distribution of 
mechanically stable configurations, and show that it is related to the existence of a peculiar infinite 
temperature critical point. 


Although systems governed by dissipative interactions 
do not obey equilibrium statistical mechanics, there have 
been several attempts to describe such systems with ef¬ 
fective equilibrium-like theories. A paradigmatic exam¬ 
ple is the problem of amorphous packings of frictional 
grains, for which an effective thermodynamic was pro¬ 
posed by Edwards and coworkers [1-5]. This approach 
relies on the basic assumption that all mechanically sta¬ 
ble packings of grains occupying the same volume have 
the same probability. This is expected if the system 
is repeatedly perturbed with “extensive operations” [3], 
like a shaking of the grains followed by a fast relaxation 
to a blocked (mechanically stable) configuration. One 
can then build an effective thermodynamics by determin¬ 
ing all mechanically stable configurations (MSCs) of the 
grains, and computing the mean values of physical ob¬ 
servables from flat averages over accessible blocked con¬ 
figurations. The predicted mean values can then be com¬ 
pared to dynamical averages obtained from a given “tap¬ 
ping” protocol which samples blocked configurations. 

Eor athermal systems in which an energy is defined, 
Edwards’ prescription can be formulated as follows [6]. 
One postulates the existence of an effective temperature 
^Ed = /^E^d the probability of a blocked config¬ 

uration C of energy E{C) takes the form 

P(C) = 2e-fedB(C) jr(c), (1) 

where Z is a generalized partition function; TiCA) = 1 if 
C is a MSG and E{C) = 0 otherwise (only blocked config¬ 
urations have a nonzero probability). This constraint is 
non-Hamiltonian, in the sense that it gives a zero prob¬ 
ability to (mechanically unstable) configurations having 
a finite energy, whereas they would have a finite prob¬ 
ability in canonical equilibrium. At first sight, Eq. (1) 
looks like a harmless generalization of equilibrium statis¬ 
tical mechanics, by simply restricting the set of accessible 
configurations. Eor instance, introducing an upper bound 
|a;i| < ATmax on harmonic oscillators Xi does not deeply 
affect their statistical properties. However, the nontrivial 
point is that the constraints arising from MSG are often 
much more complex than simply introducing a bound on 


individual variables, in particular the constraint of me¬ 
chanical stability may itself introduce strong correlations 
in the system. Notice that, beyond volume and energy, 
one may also take into account other quantities when 
building an Edwards-type thermodynamics [7-11] (e.g., 
the stress tensor). 

Several attempts have been done to test the Edwards 
scenario, not only in packings of grains, experimen¬ 
tally [12-15] and numerically [16-20], but also, in ab¬ 
stract models like spin and lattice gas [21-26], and in 
glass and spin-glass models [6, 27-30]. Typically, one 
uses a specific tapping protocol to sample blocked states, 
and compares the dynamical average of the observables 
to the thermodynamic averages obtained from Eq. (1). 
Although it has been shown explicitly in some cases that 
Edwards approach is not exact [25, 26], it is generally 
believed to be a reasonably good description in many 
cases [5]. The main difficulty with the Edwards measure 
Eq. (1) is that the partition function Z, from which all 
thermodynamic quantities can be derived, is very compli¬ 
cated to compute due to the complexity of the function 
AF{C) characterizing blocked states [31-35]. Standard ap¬ 
proaches are then either to consider abstract models [21- 
26], which are far from any realistic system but sim¬ 
ple enough to allow for an explicit solution, or to resort 
to mean-field [36] or more involved [31] approximations, 
which still capture part of the interesting phenomenol¬ 
ogy, but (at least partly) miss relevant information about 
spatial correlations in the system. 

In this Letter, we introduce a realistic model in which 
Edwards thermodynamics can be computed exactly. We 
investigate a one-dimensional model of frictional blocks 
connected by harmonic springs, subjected to a tapping 
dynamics. Due to the one-dimensional geometry, statis¬ 
tical properties can be computed semi-analytically in the 
thermodynamic limit using a transfer operator method. 
Our numerical simulation and theoretical results lead 
both to an infinite temperature critical point, with a cor¬ 
relation length diverging linearly with the stored energy 
density -a directly measurable quantities in numerical 
simulations. We analytically confirm these results using 
a Gaussian approximation for the joint probability dis- 
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tribution of spring elongations, and further characterize 
this critical point in terms of the divergence of energy 
and length fluctuations. 

Simulations- Our model is represented by a one¬ 
dimensional chain of massive blocks connected by N har¬ 
monic springs sliding on a horizontal plane [37-41]. Each 
particle is subjected to dry (Coulomb) friction. The 
position of the i**'-mass is denoted as Xi. The effect 
of dry friction is twofold: when a block is sliding it is 
subjected to a dissipative force proportional to the dy¬ 
namic friction coefficient, fi^diss = —Atdwg sgn(ii); when 
at rest, in order to make it move one has to apply a force 
larger than the static friction \fi\ > rngfis- We denote 

= Xi — Xi-i — Iq the elongation of the spring connect¬ 
ing block i and block i — 1; Iq is the rest length of the 
spring. Taking into account an external force the 
total force exerted on a block is fi = — ^i) + /“*, 

with k the spring stiffness. We write the equation of mo¬ 
tion in dimensionless form using the variables t = t/xQ, 
i = x/{gT§), = fry (mg) and Iq = lo/{gT§), with 

To = i/k/m. We have run simulations of a chain with 

-|- 1 = 256 blocks with open boundary conditions - 
see Supplementary Material (SM). We checked that, over 
the range of forces explored, no finite size effects appear 
changing the size from 64 to 2048 blocks. Dropping the 
tildes, the dimensionless equation of motion reads 

Xi = -Aid sgn(a:i) -b x^+l + Xi-i - 2xi + fr\ (2) 

while the condition to start motion becomes simply 
l^i+i — + fr*‘\ > t^s- We identify the “blocked” con¬ 

figurations as those that in absence of external force are 
mechanically stable: V i, = 0 and < fig- We 

then define the following tapping dynamics: the external 
forces /“* are switched on in Eq.(2) and act during a 
given period of time r, after which they are switched off 
and the system relaxes to a MSC. This procedure, that 
we call driving cycle, is repeated a large number of times 
to sample MSCs. At each cycle, the forces are drawn 
(randomly for each site i) from a distribution 

p(/r‘) = (1 - p) <5(/r‘) + -7= (3) 

V ZTTCr^ 

A driving protocol is determined by fixing the parameters 
p and cr. For a given protocol, one can then vary the 
intensity F and duration r of the driving. Each MSC is 
characterized by the typical value of the energy stored 
by the springs e = (1/2A^) C?- For each tapping 

protocol, the average energy e{F, t) of the MCSs is found 
to depend only on the product Ft (see SM). 

To characterize the MSC we focus on the correlation 
function C{\i—j\) = {£,i£,j) between the elongations of the 
springs at position i and j in the chain. Since this func¬ 
tion is trivial {C{\i — j\ = CoSij) in a thermal harmonic 
chain at all temperatures, any appearance of correlations 
is a signature of the unusual statistics associated with 



FIG. 1. a), b),c): correlation functions C(|i — jj) for different 
tapping protocols, while in each panel the different curves 
correspond to different energie values, a) F G [20,128], p = 
0.3, a = 0; b) E G [20,140], p = 0.8, a = 0; c) E G [20,128], 
p = 1, a — E/4, r = 60 in all simulations. Inset: C{\i — j\) 
vs. \i — j\/£{e), showing a good collapse of the different curves, 
d): correlation length £{e) as a function of the energy density 
e of the MSC, for the different tapping protocols shown in a) 
diamonds, b) circles and c) squares, showing a linear increase 
£{e) oc e with a protocol-dependent slope. 


the non-Hamiltonian constraints. Correlation functions 
C{\i — j\) measured for different tapping protocols are 
shown in Fig. 5. For a given tapping protocol we find 
that the extent of correlations increases when the average 
energy of the MSC increases. We extract the correlation 
length f(e) for each case as the distance (i.e., number of 
springs) j* — jj at which the measured correlation func¬ 
tion decays below a conventional threshold C* = 0.2. 
The insets of Fig. 5a), 5b) and 5c) show the collapse of 
the correlation function when the x axis is rescaled with 
£{e). Fig. 5d) shows, for all protocols studied, the correla¬ 
tion length growing linearly with the energy £(e) ^ e; the 
higher the energy, the more the system is correlated. This 
result may look puzzling at first sight, since commonly 
the larger is the energy the smaller is the extent of corre¬ 
lations. The key point to understand the physics of our 
system is the non-Hamiltonian nature of the constraints 
defining the MSC. If for a certain i we have ^ 1 (a sit¬ 
uation that is typical of a high energy MCS), then, to ful¬ 
fill the frictional constraint [■Ci+i — Cil < Esj f,i+i must be 
close to y. The same argument relates ^i +2 to fi+i) aud 
so on. Therefore, correlations between spring extensions 
build up in the MSC. We also compare the correlation 
functions characterizing the blocked states and those at 
the end of the driving phase (when the force is switched 
off) and show that the spring-spring correlations are not 
due to the external driving, but come entirely from the 
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FIG. 2. Distributions of springs elongations P{^) at different 
temperatures, on a semi-log scale. Points represent P{^) in 
MSC sampled via tapping, full lines represent P(^) obtained 
from Edwards theory, computed at the temperature T^d yield¬ 
ing the same energy density e as the tapping. a): uncorrelated 
regime (^(e) < 1) where e ~ Tna; Tfed = 0.0008 (blue squares) 
and 0.0026 (orange circles), b): correlated regime (£(e) > 1) 
where e ~ ^T^d', Thd = 1.19 (orange circles) and 3.84 (blue 
squares). 


constraint imposed by static friction (see SM). 

It is interesting to notice that the distribution of spring 
lengths is Gaussian at all energies. This is shown in 
Fig. 2, both in the case i{e) < 1, namely when there 
is no correlation between springs elongations, and in the 
case i{e) > 1 . 

Effective theory: transfer operators- Given that MSG 
are defined in our model by the constraint 
Vf, from Eq. (1) the probability of a configuration ^ = 
(Cl, ■ • ■ ,Cjv) reads 


N 

P(|) = e-fe- n - ^*1)- (4) 

i=l 

All the properties of the system can be obtained from 
the partition sum Z = dfi.. -P(C)- Using the 

change of variables Ci = MsCi, tire partition function de¬ 
pends only (up to an irrelevant prefactor) on the product 
/^Ed/ig! hence, all thermodynamic quantities are functions 
of TEd/Ms- We consider periodic boundary conditions 
for the chain, without imposing any constraint on its to¬ 
tal length. For convenience, we fix the rest length to 
^ 0 = 00 , allowing us to take as domain of integration 
Ci S (— 00 , 00 ), while avoiding crossings of masses. 

Using Eq. (4), we have Z = Tr(T^), with T an op¬ 
erator defined as T[/](x) = J^^dyT{x,y)f{x), being 
T{x,y) the symmetric function: 

T{x, y) = ©(^ig - lx - y\) ( 5 ) 

The operator 'T has a maximum positive eigenvalue 
Amax(/3Ed, Ms), which Can be computed numerically dis¬ 
cretizing the domain of C, and using a complete or¬ 
thonormal basis in Lf. All relevant thermodynamic ob¬ 
servables are computed in the same way (see SM). The 
free-energy is obtained as / = /^gd ln(Amax(;dEd, Ms)) 



FIG. 3. Main : Energy density e of MSG as a function of 
Tes, from transfer operators (small red diamonds) or Gaus¬ 
sian approximation (full line), and as function of dissipated 
energy Cdiss for two tapping protocols: 1) p = 0.3, a = 0 
(orange circles); 2) p = 1, cr > 0 (blue triangles). For 
Tfed <S M^ one finds an “equilibrium-like” regime, e ~ Tna; for 
Thd S> M^ the behavior is e ~ VUnd. Top inset: Gorrelation 
length from exact calculation (transfer operators): the be¬ 
havior £(TEd) ~ UTfed ~ e is clear when TEd M^- Bottom 
inset: same symbols data sets of main panel are collapsed by 
just rescaling the x-axis: up to a protocol-dependent prefactor 
we have Cdias ~ Thd. 


while the energy reads e = i9(,0Ed/)/9/3Ed = 

-(Amax|9T/9/3Ed|Amax)/Amax- M the following, we com¬ 
pare results from theory and simulations by tuning the 
temperature TEd such that the energy e takes the same 
value as in the numerics. 

The behavior of energy as a function of TEd from the 
transfer operator approach is shown in Fig. 3. We find 
two regimes separated by a crossover that depends on 
Mg: for TEd ^ Ms there is an “equilibrium-like” regime 
where e ^ TEd while for Teu Ms one finds e ^ -s/lEd- 
The transfer operator approach allows us to compute also 
the probability distribution p{f) of the elongation of a 
single spring (see SM). The theoretical result for p{f) is 
compared in Fig. 2 with the one estimated numerically 
from the MGSs, showing good agreement. We find that 
p{f) is Gaussian in all regimes, even when correlations 
are present. 

We also compute the correlation Cdz — j|) = {fifj), 
which is very close to an exponential form for all val¬ 
ues of TEd, C{\i - j\) oc e I® ilA(/3Ed,/i) (ggg g]y[)^ When 
TEd Ms both the correlation length £ and the energy e 
grow as VTeu (see inset of Fig. 3), implying £ ^ e. We 
thus recover from Edwards thermodynamics the scaling 
behavior of correlation length with energy observed in 
the simulated tapping dynamics. This is a remarkable 
success of the Edwards approach for this system. Gon- 
versely, there is almost no correlation between neighbor¬ 
ing springs {£<!), in the “equilibrium-like” low energy 
regime (e T^d)- 

We further show in the second inset of Fig. 3 that a di- 
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rect measure of T^d (within a protocol dependent factor) 
is obtained from the dissipated energy per tapping cycle 
and particle, Cdiss = Mdi/o sgn{xi{t))xi{t)dt). Indeed, in 
the simulations ejiss is found to have the same scaling 
with e as the temperature obtained within the transfer 
operator approach (i.e., e ~ Cdiss if e ^ e ~ -y/cdiss if 
e ^ /Tg). The dissipated energy can therefore be inter¬ 
preted as the analog of the thermal energy that allows 
the system to sample the configuration space [42]. 

Effective theory: Gaussian ansatz- To obtain ap¬ 
proximate analytical expressions for the thermodynamic 
quantities, we replace the Heaviside function in Eq. (4) 
by a Gaussian function, 

0 (Ms-|6-ri-6l) ^ ^ exp (6) 

yielding 2 oc f with an effective Hamiltonian 


S{0 


1 

2 


N 


PEd ^ 




(7) 


This effective Hamiltonian corresponds to a positive def¬ 
inite quadratic form S{^) = A^/2, where H is a sym¬ 

metric real Toeplitz matrix (see SM). The matrix A can 
be exactly diagonalized, yielding analytical expressions 
for energy, entropy, correlation function and correlation 
length. The mean energy per particle reads 


e(rEd,/rs) 


1 Ms^Ed 

2 ^2rEd + til ’ 


( 8 ) 


from which we recover that the crossover point between 
the behaviors e ~ TEd and e ~ \/T-^d is T^d ~ Ms- 
Fig. 3, Eq. (9) is compared with the result from the 
transfer operator Eq. (5), showing a semi-quantitative 
agreement. We also find that the correlation function 
is ~ with £(TEd,Ms) such that 

i ~ a/T/Ms in the limit Tsd Mg. This result can be 
recovered from a field-theoretic viewpoint, by taking a 
continuous limit in Eq. (7), yielding 


-S'K] oc 


dx 


1 

2 





(9) 


with a “mass” term nff = 2ff^/3-Ed- The correlation func¬ 
tion of such a Gaussian field theory is known [43] to 
be {^{x)^{y)) ~ so that we recover a corre- 

lation length £ ~ . This field-theoretic formulation 

confirms the presence of an infinite temperature critical 
point, since the mass term goes to zero at infinite tem¬ 
perature. To inspect the critical exponents associated 
with this critical point, we study the fluctuations of the 
total energy of the chain, 6E = ^ ~ (CD) and the 

fluctuations of its total length S£ = We find (see 


SM) that the variance of both energy and length diverge 
linearly with temperature (or, equivalently, as f^), 


N 


N 


Tfed • 


( 10 ) 


Finally, we compute the entropy density s = —df/dTEd- 
We find that it saturates at high temperature to a fi¬ 
nite value, limTEd->oo s(TEd, Ms) = 5 In 2 -I- In /Xg . This 
saturation results from the presence of long-range cor¬ 
relation at infinite temperature. This can be confirmed 
by contrast, computing the “mean-field” entropy density 
Smf = —f d^p(^) Inp(^), with p(^) the distribution of 
a single spring elongation We find that Smf) which 
discards correlations, diverges like InTsd at infinite tem¬ 
perature (see SM), at odds with the saturation of the 
entropy s. 

Conclusions- The present study provides a clearcut 
example of how an effective thermodynamic theory can 
successfully describe an athermal dissipative system. 
The most remarkable difference between standard equi¬ 
librium thermodynamics and the effective theory we have 
presented is the presence of an infinite temperature crit¬ 
ical point, with an associated divergence of the corre¬ 
lation length as £ ~ (or I ^ e). As seen in the 

field-theoretic formulation Eq. (9), this infinite temper¬ 
ature critical point results from the long-range correla¬ 
tion generated by static friction in the blocked states. 
The difference with standard equilibrium systems is that 
the gradient term in the effective Hamiltonian does not 
come from an energetic interaction, but from a non- 
Hamiltonian constraint. Its coefficient is strictly tem¬ 
perature independent, while the coefficient of the en¬ 
ergetic term scales inversely with temperature. While 
temperature-independent terms could also be present at 
equilibrium (e.g., entropic constraints such as excluded 
volume), they are usually purely local and do not in¬ 
volve gradient terms. Hence, in spite of its simplicity, 
our model exhibits a phenomenology clearly distinct from 
that of equilibrium systems, and the field-theoretic for¬ 
mulation suggests that the results should be quite ro¬ 
bust to changes in the details of the model. Future work 
should investigate this issue in more details. 
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Numerical simulations 

Details on the simulations 

The tapping protocol is defined in Eq. (3) in the main 
text. The motivation for introducing (annealed) disor¬ 
der in the tapping protocol is that, when taking p = 1 
and cr = 0 (i.e., without any disorder in the driving), 
we found undamped waves which travel across the chain 
causing an undesired and artificial collapse of the “ensem¬ 
ble” of visited MSC to a very small and specific subset. 
As mentioned in the main text, for the numerical simu¬ 
lations of the tapping dynamics we consider open bound¬ 
ary conditions. This choice is done in order to avoid 
constraints on springs which may induce trivial correla¬ 
tions. For instance if one fixes then also 

has = (L^) = which is a constraint on 

the two point correlation function. With open boundary 
conditions the first and the last blocks are connected to 
a single spring and their equations of motions are respec¬ 
tively 

xi = -pdSgnii -f {X2 - xi- lo) + /r* 

Xn+ 1 = -lJ.dSgD.XN+1 - (XN+I -XN - lo) + /w+1 ( 1 ) 


Note that in order to have N springs with open boundary 
conditions we need N + 1 blocks. 


Dependence of the energy on driving parameters 

The average energy of MSC depends in principle on the 
four parameters which charaterize the driving cycle: the 
duration r of forcing, the intensity of the mean applied 
force F, its standard deviation a as well as the fraction 
p of pulled particles. The parameters p and a define the 
tapping protocol, while F and r characterize the intensity 
of the driving. We have found that for a fixed tapping 
protocol, the energy e of the dynamically reached MSC 
does not depends separately on F and r, but is uniquely 
determined by the product Ft. This is clearly shown by 
the set of data reported in Fig. 4. 


One-cycle experiment: test for the origin of friction 
induced correlations 

In the manuscript we presented data on the “spring¬ 
spring” correlation function obtained in the me¬ 

chanically stable configurations (MSCs). Due to the non- 
random type of the external force F used to sample dif¬ 
ferent MSCs the reader is allowed to wonder if the cor¬ 
relations we measure in the MSCs are really induced by 


u 

C 0.2 ' 


2000 4000 6000 8000 10000 12000 14000 

Ft 


FIG. 4. Energy per spring stored in mechanically stable con¬ 
figurations (MSC) as function of Er with F the pulling force 
and r the duration of the pulling. Full circles: F is increased 
at fixed r; Empty squares: r is increased at fixed F. 


friction at 100%, or a part of them is induced by the same 
F. The numerical evidence that it is not the case comes 
from the “one-cycle driving experiment” we are going to 
discuss. The one-cycle driving experiment is analogous 
to the tapping dynamics discussed in the main text apart 
from one thing: at the beginning of each driving cycle the 
positions of blocks are randomized. More precisely all the 
Xi are reset to Xi = i Iq + dx, where dx is normally dis¬ 
tributed and Iq is the rest length of springs. During the 
elementary driving cycle of the dynamics a fraction p of 
the particles is pulled for a durantion r with a constant 
force F, the force is then switched off and the system 
relaxes to a MSC. We average over many iterations two 
correlation functions: a) the forced correlation 
which is measured just before switching-off the force at 
r; b) The correlation in the mechanically stable config¬ 
urations reached after the quench, (^i^j)MSC- Our aim 
in studying {£,i£,j)F is to see how much of the correlation 
is provided by the external force. On the other hand we 
must also be sure that when measuring {f,if,j)F there is 
no memory of the correlation in (CiCi)MSC at the previous 
step. It is for this reason that at the beginning of each 
driving cycle the positions of block are randomized. In 
this way we make ourselves sure that if some correlation 
is measured in (^i^j)_F, it comes solely from the external 
force. The result is illustrated in Fig. 5. The continuous 
line is {£,i£,j)MSC while the (red) circles represent {^i^j)F' 
being the latter identically zero we can conclude that no 
correlation is induced by the force and that the corre¬ 
lations we observe in the system come solely from the 
constraint imposed by static friction on the MSCs. 
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FIG. 5. Spring-spring correlation function (red cir¬ 
cles) in presence of the external force and (contin- 

uons line) measured in the mechanically stable configurations 
reached after qnench. Data are obtained within a “one-cycle- 
driving-experiment”: an external constant force F = 100 is 
applied for a duration r = 100 pulling a fraction p = 0.25 of 
the blocks. The number of blocks is N = 256. At the be¬ 
ginning of each driving cycle the positions of the blocks are 
randomized. 



0 5 10 15 20 25 30 35 

i-j 

FIG. 6. Main: spring-spring correlation function C{i — j), 
where i and j labels different springs {i > j), computed ac¬ 
cording to Eq. (8) at different temperatures. From left to 
right: T = 0.185,0.341,0.631,1.166,2.154,3.981,7.356. At 
odds with what usually happens the extent of correlations in¬ 
creases as the temperature is increased. Inset: C{i — j) vs 
i — j in semi-log scale: all correlations decay exponentially. 


Transfer operators 

The core of the exact solution of the effective ther¬ 
modynamics of our model is represented by the spectral 
properties of the linear operator T, acting in the space 
of square integrable functions T : t L^. The 


operator T is defined as follows: 

/ OO 

dy T{x,y)f{y), (2) 

-OO 

with 

T{x,y) = \y - x\)e~^'^\ (3) 

As can be easily checked / dxdy\T{x,y)\^ < oo, so that 
T is a Hilbert-Schmidt operator, and hence compact and 
bounded, which guarantees the existence of a maximum 
positive eigenvalue. The kernel is symmetric and real, 
T{x,y) = T*{y,x), so that the operator is also self- 
adjoint, which guarantees the spectral theorem to hold: 
a complete set of eigenvalues and eigenvectors exist for 
our operator [ 1 ]. 

Correlation function 

By means of the transfer operator approach we can 
compute several quantities of interest: the one that is 
most relevant for our study is the correlation between 
the elongation of far apart springs. Let us consider for 
instance the elongations and with j > i. In terms 
of our effective theory the correlation between these two 
elongations can be written as: 

(c.c,) = 1 / da...eivT(ei, 6 )---ne.-i,coe.ne.,e.+i) 

(4) 

which in terms of operators becomes 

iUj) = , ( 5 ) 

where V is the diagonal operator defined as: 

V[f]{x) = xf{x). (6) 

Exploiting the fact that the operator T has a complete 
spectrum of orthonormal eigenvectors we can use the 
completeness relation 1 = X]al®)(®li where |a) is the 
eigenvector of T relative to the eigenvalue a, to sim¬ 
plify Eq. (5). One can write the trace in Eq. (5) as 
Tr(T) = then insert between each pair of 

neighboring operators the completeness 1 = 

The result is 

( 60 ) = 

^ {a\T^-^\b){b\D\c){c\T^-^\d) 
{d\D\e){e\T^-^+^\a) 

= (7) 

a,b 














which in the thermodynamic limit becomes 


lim 

N—¥<yD 




\N 
^max 




E 

b 

E 




|(&|i?|A„,ax)|' 


Ab 




/ oo 

dx xf,ix)f,^^Jx) 

-OO 


( 8 ) 


where f^{x) is the eigenfunction related to the eigenvalue 
b. The correlation functions obtained at different values 
of the Edwards temperature T^d by means of a numerical 
determination of eigenvalues and eigenvectors of T can 
be found in Fig. 6 : at all temperatures, we numerically 
observe an exponential decay, meaning that the second 
largest eigenvalue dominates the sum in Eq. ( 8 ). 


Energy 


The average energy of the harmonic chain with N par¬ 
ticles is defined as: 


eAr(/3Ed) 


dP'Ed 


(/^Ed/) = - 
1 


1 d 


HZ) 


N dp-Ed 

^ Tr[r'^(/3Ed)] 


fV Tr[r^(/3Ed)] dpdid 

1 8T 


(9) 


where in the last line of Eq. (9) we passed the deriva¬ 
tive sign d/djd'Ed under the trace sign. The “matrix el¬ 
ements” of the operator dT /9/3Ed are simply defined as 
the derivative with respect to /3Ed of the matrix elements 
of T: 


dT 

dfi-Ed 


{x,y) 


x'^ + 
4 


exp 



( 10 ) 


Exploiting again the definition of the trace Tr[T^] = 
^^(a|T^|a) and by inserting the completeness relation 
1 = |a)(a|, we can write 


eAr(/?Ed) 


1 




Tr[r^(/3Ed)] 
-^^Tr[r-^(^Ed).l.|^] 




EaA^ 


ab 


d/^Ed 


Ea, 




df^Ed 
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which in the thermodynamic limit yields 


lim Bn iPEd) 
N—¥00 


-^(X |EL|a \ 

X \^max| p.n l^max/ 

^max C'PEd 


— / dxf,^^Jx)dp^^T{x,y)f,^^Jy), 

■max J —oo 

( 12 ) 


where f\-^Hy) is the eigenfunction related to the maxi¬ 
mum eigenvalue. 


Marginal distributions 


By means of the transfer matrix technique is easy to 
obtain marginal distributions. For instance we are inter¬ 
ested in the probability distribution of the springs elon¬ 
gation, which reads 


PnH 


lim 
N —¥00 


PnH 


1 1 

“ ab 
2^a '^a ^ 

l/w.(e)l^ (13) 


where, again, (a;) is the eigenfunction related to 

the maximum eigenvalue. This distribution, which has 
a Gaussian shape, is plotted for different values of TEd in 
Fig.2 in the Letter. 


Gaussian approximation 

An explicit computation of how the internal energy 
depends on the temperature is possible if one assumes a 
Gaussian smoothening of the Heaviside function enforc¬ 
ing the static friction constraint: 

e(Ms - |g.+i - ed) ^ ^ exp (^- (14) 

With this approximation the whole partition function 
reads 


Z = 7r-^/2 J 

= J d^i... d^N 


exp 





( 11 ) 


(15) 
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where ^ represents the vector ^ = (^i,..., and A is 
a, X X N symmetric real Toeplitz matrix of the kind 


with 


/ Ao Ai 0 0 0 Ai \ 

Ai Ao Ai 0 0 0 

0 Ai Aq Ai 0 0 

0 0 Ai Aq Ai 0 

0 0 0 Ai Aq Ai 
V Ai 0 0 0 Ai Aq J 


Aq = /3-Ed H-y 

K 



(16) 


Gaussian integral formula, and all eigenvalues of A con¬ 
tribute to the free-energy (while only the largest eigen¬ 
value of the transfer operator contributes to the free en¬ 
ergy in the thermodynamic limit). 

The energy per degree of freedom can also be com¬ 
puted: 

_ d{P'Edf) _ _1_ 

a/^Ed 2A ^ ^ [l _ cos (^)] 

1 _dfc_ 

2 Jq /3Ed + ^ [1 - cos(27rx)] 

_ Ms _ 1 Ms^Ed 

" Vp Ed + /3Vs' 2 V2rEd + M^’ 

( 21 ) 


For this kind of matrix the eigenvalues can be explicilty 
written as 


Afc(/3Ed, Ms) = Aq -I- 2Ai cos 


= /^Ed H—y , 

K L 


/ 2Trk 

(2TTk\ 


(17) 


Free energy, energy and entropy 


Combining Eqs. (20) and (21), the entropy per degree of 
freedom can be simply obtained as s(/3Ed) = /3Ede(/3Ed) — 
/3Ed/(/3Ed)- One can check that s(/3Ed) —t y ln2 -|- In^s 
when Tejj —> oo. Also, energy fluctuations are given by 

{[SE]^) l + /3EdMs (QQ-. 

N 2 2/3Ed + /3idMs"’ ^ ^ 


and at high temperature they diverge as {[SE]‘^)/N ~ 
TEd, consistently with Eq. (10) in the main text. 


The free-energy reads then as 


Correlation function 


/ = - 


1 


PEdN 

In 2 

+ 


In 


1 


2/3Ed 2.j5^QN 

In 2 1 

+ 


J d^N exp 
In det(A) 




2/3Ed 2/?EdA^ 


^ln[Afc(/3Ed,Ms)]- 


(18) 


In the continuous limit 


1 ^ 

N ^ 


dk 


k=l 


(19) 


the free-energy can be exactly computed and reads: 


/ = 

2/3Ed 2/?Ed 

In 2 1 

+ 


f dk In (p-Ed + ^ [1 - cos (27rA:)]^ 
do \ Ms / 


2/3Ed 2/3Ed 


In 


/^Ed 




/^EdMs 


( 20 ) 


When the probability of a configuration of the system 
is a multivariate gaussian of the kind 


P(^) (X exp 



(23) 


we know that the value of the average correlation between 
and is 

= (AI-'),,. . (24) 

For the problem we are studying A is a tridiagonal 
Toeplitz matrix, for which the explicit formulae for the 
elements of the inverse is [1]: 




(_1)7+J 


A\-^ I/,-i(d)I/„-,-(d) 

|Aip-*+i Un{d) 


(25) 


where d = Ao/(2|Ai|), and Un{x) is a Chebyshev poly¬ 
nomial of the second kind of degree n in the variable x. 
For such polynomials, a useful identity is known: 


Note that this method to compute the free-energy is dif¬ 
ferent from the transfer operator method. Here we do 
not write the partition function as the trace of the A**' 
power of a transfer operator, but we diagonalize exactly 
the quadratic form appearing in Eq. (15). The partition 
function is then exactly computed for all N using the 


Unix) 


sin0(n -I- 1) 
sind 


where 9 is the complex number such that x 
0 < 5ft(0) < TT. In our case we have 


(26) 
cos(d) and 


(i = Ao/(2|Ai|) = I + ^2/3Ed>l, (27) 
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where the last inequality on the right is true for every 
finite value of the temperature and of the friction coef¬ 
ficient. Therefore the number 9 must be a purely imag¬ 
inary one, i.e., 9 = ia, with a > 0. For our problem 
we have therefore that the Chebyshev polynomials are 
defined according to the two following identities: 


[/„(! -I- /ig/^Ed) = 


sinh Q;(n + 1) 
sinha 


1 -I- fJ,sf3Ed = cosh(Q;) 


(28) 


Because our problem is in the continuum we are inter¬ 
ested in the limit n —> oo of the expression in Eq. (25). 
Let us notice the following: 


which scales for Tsd —>■ oo as £{^s,T-Ed) ~ V^Ed/Ms- 
The variance of the fluctuations of the total length C = 
Ef=i read 


with C{r) defined by C{\i — j\) 
we get 


N / dr C{r) 


(33) 


= i^iQ- Using Eq. (31), 




/ oo 

dr exp{—r/£) ~ Ne£ ~ NTe^ (34) 


Jo 

thus recovering the second equality in Eq. (10) of the 
main text. 


lim 

n—>-oo 


Un-jid) 

Un{d) 


C/._i(d) 


lim 

i<j 


U,.lid)Un-jid) 

Un{d) 


sinh[(n -3 + l)a] - 

n->oo sinh[(n -f l)a] 
sinh(tQ;) 
sinh(Q;) 

_ p-a(i+j) 

lim ^-—fx- 

i-^-oo 2sinh(a) 

i<j ^ ^ 


2 sinh(£-i)’ 


(29) 


where the correlation length £ in the last row is defined 
a.s a = £~^. In the limit where UsPeo ^ Ij namely 
when the temperature is large compared to the friction 
coefficient, the correlation length £ must be large and we 
can estimate its asymptotic scaling: 


cosh(£ ^) 

Ms/^Ed ^ 1 
Ms/^Ed ^ 1 


= Ms/^Ed + 1 
Ms/^Ed + 1 



Ms 



(30) 


By recalling that Ai = —l/{2fd^), we can finally write 
the two-point correlation function as: 


Spring length distribution 


The probability P{^) is a multivariate Gaussian, so 
that the marginal distribution p(^) = f d ^2 ■ ■ ■ d^NP{^) 
is also a Gaussian: 


PiO 


1 

- -- g 2cr^(TEj,Ais) 

fr(TEd,Ms)v^ 


(35) 


with 

a2(TEd,M.) = {A-%, = (36) 

V 2TEd + Ms 

The “mean-field” entropy Smf introduced in the text can 
be computed as 

Smf = - y d^ p{^) ln[p(5)] 

= ^ -I-ln[cr(rEd,Ms)VE-], (37) 

which yields the asymptotic result 

Smf(7Ed) ~ InTEd , Tfed —>■ OO . (38) 




(60) 


_ ds _ -a(j-i) 




sinh(a) 


OMs + 2 rEdMs 
2 e(Ms,TEd) 


(31) 


The fact that the entropy of Eq. 37 in the infinite tem¬ 
perature limit diverges, and is then wrong compared to 
Eq.lO in the Letter, is consistent with the fact that for¬ 
mula in Eq. 37 is approximation which does not take into 
account the correlations. 


where the correlation length is 


£{Ps, T^d) aj.(;osh(l -I- Ms/^Ed) 


(32) 
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